Screening of magnetic fields by charged Bose condensate. 
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Abstract 

The space-space component of the photon polarization operator is calculated in zero fre- 
' quency limit for a medium with Bose-Einstein condensate (BEC) of electrically charged par- 

ticles. It is found that the polarization operator tends to a finite value at vanishing photon 
3-momentum, as it happens in superconducting media. It means that magnetic fields are ex- 
ponentially screened in such a medium analogously to the Debye screening of electric charges. 
' At non-zero temperature the screened magnetic field oscillates and contains a contribution 

. which drops only as a power of distance. This phenomenon is unknown for superconductors, 

even in BEC phase and can be potentially observable. 

As is well known, electrically charged impurities in plasma are screened according to the 
Debye law [ij: 

^ q exp(-mj,r) 

47r r ' 

(N 

• where q is the charge of a test particle and rriD is the screening (Debye) mass, which is expressed 

. through the density of charge carriers in plasma. On the other hand, it is known that magnetic 

I fields are not screened in plasma because of the absence of magnetic monopoles. The last state- 

' ment has been proven in standard QED [3j and in scalar QED (SQED), when bosons are not 

. condensed. We show below that, in the presence of Bose-Einstein condensate (BEC) of charged 

' particles, magnetic fields in plasma are exponentially screened. Such screening is a manifestation 

of the well-known Meissner effect in superconductors. However, at non-zero temperature new 
phenomena of an oscillating exponential screening and a power law one are found here. The 
screening mass is of course different from the usual Debye one, which, in the relativistic case, is 
proportional to the temperature or the chemical potential. The charged condensate contributes 
to the photon polarization tensor by additive terms, whose amplitude is proportional to the am- 
plitude of the condensate and is the same for electric and magnetic photons. The calculations 
are done in the framework of SQED with charged scalar boson condensate but it is evident that 
similar screening exists in the presence of charged vector condensate. Moreover, the result seems 
to be applicable to non-Abelian gauge theory as well. These problems will be studied elsewhere. 

The electric screening in plasma in the presence of charged Bose condensate was first studied 
in refs. [3] and recently in relativistic theory in refs. [3 [6l [71 [8] j see also [9] and references therein. 
It has been discovered that the photon polarization operator in plasma is singular at zero photon 
3-momentum, having poles proportional to iSiiT] and 1/k [6l|8], which lead to an oscillating 
and power law screening behavior. Such screening behavior is in good agreement with the earlier 
papers [Ij. Magnetic, as well as electric, screening with charged BEC at zero temperature was 
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considered also in ref. [TU], where the superconducting behavior as a response to electromagnetic 
field was observed. 

The Fourier transform of the Maxwell equations in plasma is modified as: 

[kPkpg^"' - k^^k" + W{k)] A^{k) = J^'ik) , (2) 

where Af^ is the electromagnetic potential, k^ is the photon 4-momentum, J'^{k) is the 4- vector 
of an electric current, and H^^ is the polarization operator of photon in medium. The latter can 
be expressed through the distribution functions of charged particles and antiparticles in medium, 
f{q) and f{q), where q is the particle 3-momentum (see below). 

The expression for the polarization operator can be found in many textbooks on field theory at 
non-zero temperature. For details of calculations we address to our work [6j, where the notation 
is the same as here, to refs. [8| and [2] and references therein. In the lowest order in electric 
coupling, the photon polarization operator induced by interactions with charged bosons has the 
form: 



(27r)3S 



Ub + Ib) 



1 {2q 
2 



k)^{2q-k)^ , 1 {2q + k)^{2q + k), 



{q - kf 



m 
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{q + kf 



m 



B 
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g2 + m o and we omitted for brevity the arguments of / 



where index B denotes bosons, E 
and /. 

In thermal equilibrium the distribution functions of bosons and anti-bosons, respectively /b 
and fs, take the form: 

AC) _ 1 
Jb 



exp[(^ -mB)/T\-l 
1 



B 



(4) 
(5) 



exp[(^ + mij)/r] - 1' 

where the chemical potential of bosons is equal to their mass: fi = uib and = —ruB for anti- 
bosons. The constant parameter C describes the amplitude of the condensate. Note that bosons 
condense when their chemical potential reaches the maximum allowed value equal to the boson 
mass. 

The time-time component of the polarization operator, IIoo, was analyzed in detail in refs. 
[6l[8]. Such a quantity, in the limit of vanishing photon energy, a; = 0, and small magnitude of 
the photon three-momentum, — > 0, is relevant to the screening of electric charge. The screening 
of magnetic fields is connected with the space-space components, Iljj, which, in the homogeneous 
and isotropic case, can be presented as: 



a{k) 6ij 



+ h{k)'^' ' 



(6) 



To calculate a{k) we first multiply H^j by 6ij and, second, by kikj. Correspondingly we obtain: 



2a{k)+b{k) 

m 

where kq = kq cos 6. Hence 
aik) = 



167r3 
0, 



2kq - A;2 2kq + k^ 
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If only the condensate term (delta- function) is retained in distribution function ([H, the inte- 
gration over d^q is trivial and we obtain: 



2 2 

e m,j , 



(10) 



Since a{k) does not vanish at A; = but tends to a constant, the asymptotic behavior of magnetic 
fields at large distances from the source would be given by the exponentially decreasing Yukawa 
potential, i.e. essentially the same as the Debye screening. The existence of the magnetic mass 
in a medium with a BEC was found in ref. [7j, but the screening effects were not emphasized 
there. In the subsequent paper by the same authors [11] it was shown that below certain critical 
value magnetic field is expelled from the condensate, while for larger values it penetrates the 
condensate within flux-tubes. 

On the other hand, if we consider the standard Bose distribution with ^ < niB, when no 
condensate is formed, the distribution has the usual form: 



1 



B 



:ii) 



then Ilij{k) vanishes as k in the limit k 
easily taken and we find: 



eMiE - li)/T] - 1 ' 

0, as expected. Indeed, the integrals over angles is 



a{k) 



167r^ 



q^dq 
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4q2 - k^ 
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2q + k 
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Expanding the integrand in in powers of (k/q) we find: 



a{k) 



247r2 



(13) 



As expected, a{k) ^ k"^ and magnetic fields are not screened. 

Note that expression (jl3p is singular in the limit niB = 0, since the integral diverges as 
1/q'^ at the lower bound of integration, q = 0. Even if there are no massless charged particles, 
somewhat similar singularity exists for massive ones if /i = and /c — )• 0, as it was found in 
ref. [6]. Such singularity leads to a term ~ 1/k in the time-time component of the polarization 
tensor. Similarly, the same singularity leads to a non-standard contribution to the space-space 
component: a{k) ~ k, when A; ^ 0. This can be demonstrated by evaluating the integral in (|12|) 
with the distribution function (jlip where fi = ttib- The additional, harder in k, terms come from 
the integration region where q is small, i.e. comparable to k. It is convenient to separate the 
integral into two parts < q < k/2 and k/2 < q < oo. In the first part we introduce the new 
integration variable x = 2q/k, so < x < 1. In the limit of small k the energy can be expanded 
as Eb ^ tub + k'^x^ /SniB- At small q the distribution function is infrared singular: 



exp 



Eb - ruB 
T 



1 



2mBT SuibT 
g2 k^x^ 



(14) 



Usually this singularity is not dangerous because it is canceled by the integration measure, ~ g^. 
However, the logarithmic term behaves as k/q for q > k and as q/k for q < k. Thus the integral 
is finite, but it does not vanish as /c^ when k ^ 0. 
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The first part of the integral with q < k/2 can be taken analytically and we obtain: 



e^kT 
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There is also another contribution coming from the part of the integral with q > k/2. As k ^ 0, 
the second part of the integral, k/2 < q < oo, gives: 



(«)^^.^ ^^^^ 
«2 (k) = 
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such that the total contribution is: 



a(^\k) = a['\k) + ai'\k) = ^k. 



(16) 



(17) 



For small k this term could dominate over the usual /c^ term and would change the screening 
behavior. 

In the transverse gauge, k^A^^ = —kjAj = 0, equation ([2]) takes the form: 

'k^ + a{k)]A, = J,{k), a{k) = a^'\k) + a^^\k) (18) 
whose solution in coordinate space can be written in terms of the Green's function as: 

A{x) = J d^yG{x-y)J,{y)- (19) 
The asymptotic behavior of G is determined by the integral (see e.g. ref. [8]): 

roo dkk {e^^^' — e~^^^^ 



G{r) 



47r2 



r Jo 



fc2 + a{k) 



(20) 



According to the result obtained above, a{k) may contain odd terms in k and hence the integral 
along the half real fc-axis cannot be extended to the whole real axis, so we present the denominator 
as half of sum and difference of even and odd function as following: 



fik) = [fik) + f{-k)]/2 + [f{k) - f{-k)]/2 . 
Since a{k) = k"^ + e^m^ + e^Tk/lQ, eq. (j20p can be rewritten as: 



(21) 
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dkk (e''"'' - e"*'''') {k^ + e^ml - e^Tk/lQ) 
(fc2 + e2m2)2_e4r2A;2/256 



(22) 



The integral of the even part may be transformed, as usually, into the integral along the whole 
real axis and after closing the contour in the upper (for e*'^^' ) or lower (for e~'^^^' ) half-plane we 
express the result through the residues in the corresponding poles in the complex fc-plane at: 



fc -±t^emc ^Q24^ 32 



(23) 
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If 771(7 > e^T /2>2, the resulting screened potential would be exponentially cut with superimposed 
oscillations. For e^T <^ 32 mc, the Green function takes the form: 

G{r) ~ exp(— emc?') cos(e^rT/32) . (24) 

In this case the spatial damping scale is much shorter than the oscillation scale. However, if 
eT ~ mc the scales are comparable. These oscillations remind the Friedel oscillations [12j in 
fermionic plasma, but physics is different. 

The contour of the integration of the second, odd in k, part can be closed in the upper or lower 
quadrant of the complex /c-plane. So in addition to the poles in these quadrants we need to include 
the contributions from the integrals over the imaginary axis. The latter produces a power law 
asymptotics of G{r). If the condensate is absent, but /i = ms, then at small k: a{k) ^ e^/cT/16, 
see eq. ()17p . Correspondingly at large r the Green's function drops as: G(r) ~ 8/(7r^e^r^r). This 
asymptotic behavior is realized when r > 1/T. In the presence of the condensate the Green's 
function acquires an additional constant term e^m^ (llOj) . In this case the contribution of the 
integral over the imaginary axis of k gives G ~ T/(16 e^vr^r^m^). 

In ref. |13j it was found that W bosons, which may condense in the early universe in the 
presence of a large lepton asymmetry would form a ferromagnetic state with spins aligned 
on macroscopic scales. As it was pointed out in the quoted paper, the ferromagnetic alignment 
could be destroyed by screening of magnetic interactions. By comparing the screening length 
A = 1/emc with the typical inter-particle distance d = C~^^^, one can roughly check whether 
screening is effective or not. It follows that screening is effective unless vtib/C^^^ <^ 10~^. Thus 
it seems hardly possible that the ferromagnetic state could be realized in the broken phase of the 
electro- weak sector, where the boson mass, m^, is essentially determined by the Higgs VEV. On 
the other hand, screening may be ineffective in the unbroken phase, where is much smaller 
and determined by radiative corrections. 

The exponential screening of magnetic fields and their oscillations in massless scalar electro- 
dynamics were mentioned in ref. [H]. They were justly prescribed to an inaccurate treatment of 
the perturbation theory with massless charged bosons. If one takes into account the temperature 
corrections to the charged boson mass by the resummation of the perturbation series, the effect 
disappears and the polarization operator vanishes as k'^ at small k in accordance with the expec- 
tations. However, in our case the temperature corrections to the mass do not help to avoid the 
effects considered here. According to kinetic treatment, the condensate appears when the boson 
chemical potential, /i, is equal to the boson mass. So the temperature corrections to the mass are 
canceled by equal correction to the maximum value of /j, and the singularity is not eliminated. 

The exponential screening in the considered system indicates that the system is similar to the 
superconductor in strongly coupled BEG phase. However, the new effects of power law decrease 
of the magnetic field and its oscillations are new and quite surprising. Possibly the same or 
analogous phenomena exist in the usual superconductors and may be observable. 

The phenomenon of the power law magnetic field screening may have interesting applications 
to compact stars with magnetic fields. This effect could also give rise to an infrared cutoff in 
non-Abelian gauge theories, which suffer from very strong infrared singularities [16j. Moreover, 
the infrared divergences in non-Abelian theories could enforce or facilitate the formation of Bose- 
Einstein condensation of charged bosons. 
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